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The  work  presented  in  this  report  was  performed  at  the  General 
Electric  Company,  Space  Division,  Space  Sciences  Laboratory,  King  of  Prussia, 
Pennsylvania.  The  work  was  done  under  the  auspices  of  the  Office  of  Naval 
Research  with  Mr.  J.  A.  Satkowski  at  technical  monitor. 

Dr.  Charles  H.  Marston  was  responsible  for  channel  design  and 
theoretical  investigation.  Mr.  E.  Tate  was  responsible  for  experimental 
operation  and  diagnostic  evaluation.  Messrs.  W.  Frey,  F.  McMenamin  and  G.  Fecik 
provided  technical  assistance  in  facility  operation.  Mr.  V.  Tilli  was  responsible 
for  mechanical  design  and  detailing  of  the  new  MHD  channel.  Mr.  Tate  and  Mr. 

L.  D.  DeDominicis  were  technical  liason  with  NASA  Ames.  Dr.  Bert  Zauderer, 

Manager  of  MHD  programs,  was  the  principal  investigator. 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


I 

I 


I 


I 

I 

I 

I 

I 


TABLE  OF  CONTENTS 


Page 


SUMMARY  OF  KEY  RESULTS iv 

1.  INTRODUCTION 1 

2.  THEORETICAL  DISCUSSION  AND  DESIGN  CONSIDERATIONS A 

2.1  Introduction A 

2.2  MHD  Generator  Model 5 

2.3  Electromagnet  Coil  Model 6 

2.  A Magnetic  Field  Model 7 

2.5  Circuit  Model  for  Combined  Coil  and  Permanent  Magnet...  8 

2.6  PS  EMC  Circuit 9 

2.7  Design  Considerations 12 

3.  EXPERIMENTS  AND  DESCRIPTION  OF  HARDWARE 1A 

3.1  ONR/GE  Shock  Tunnel 1A 

3.2  GE  I MHD  Channel 15 

3.3  NASA  Ames  "EAST"  Facility 17 

3.  A Ames  I Channel  Configuration  and  Experiments 18 

3.5  Ames  II  Channel  Configuration  and  Experiments 20 

A.  CURRENT  STATUS 22 

5.  REFERENCES 23 

6.  FIGURES  AND  TABLES 2A 

APPENDIX  A A-l 

APPENDIX  B B-l 

APPENDIX  C C-l 


iii 


■ . 


SUMMARY  OF  KEY  RESULTS 


1.  Feasibility  of  self-excited  MHD  operation  has  been  shown  with  the  GE  I 
MHD  channel.  8%  magnetic  field  augmentation  using  an  initial  field  of 
0.5  Tesla  was  obtained. 


2.  The  NASA  Ames  "EAST"  facility  capabilities  as  a high  temperature,  high 
pressure  plasma  source  have  been  mapped  in  detail. 

3.  A channel  (Ames  I)  has  been  used  on  the  EAST  facility  to  check  test-time, 
aerodynamic  performance  and  construction  materials. 


4.  Design  and  construction  of  a preliminary  (Ames  II)  pulsed,  self-excited, 
MHD  generator  concept  was  completed. 


5.  Initial  experiments  were  done  on  the  EAST  facility  using  the  Ames  II 
channel  and  some  problem  areas  identified. 


1.  INTRODUCTION 


The  long  term  goal  of  this  program  Is  the  direct  conversion  of  energy 
stored  in  a solid  explosive  to  electrical  power  utilizing  an  MHD  generator. 

The  power  density  of  an  MHD  generator  is  proportional  to  the  electrical 
conductivity  and  a solid  explosive  can  produce  very  high  (>  10,000°K)  gas 
temperatures  and  hence  high  electrical  conductivity.  The  work  proposed  is 
primarily  concerned  with  the  conversion  of  energy  in  the  conducting  gas  to 
electrical  energy  and,  while  several  methods  for  explosive  heating  are  possible, 
they  are  not  considered  here. 

To  be  of  interest  as  a practical  device,  an  explosive-MHD  power  source 
should  be  compact,  self-contained  (i.e.  no  external  ausiliary  power  and  little 
or  no  preparation  for  use  required)  and  it  should  have  a long  storage  life. 
Unfortunately,  none  of  the  previous  experiments  met  these  requirements  because 
of  the  need  for  providing  a magnetic  field  for  the  MHD  generator.  To  achieve 
practical  overall  conversion  efficiencies  (ranging  from  1 to  5%)  it  is  necessary 
for  the  electromagnetic  forces  in  the  MHD  generator  to  be  of  the  order  of  the 
gas  dynamic  forces  in  the  explosively  produced  plasma.  The  latter  has  a pressure 
of  hundreds  of  atmospheres  and  for  these  conditions  the  explosive  MHD  generator 
needs  to  operate  at  magnetic  fields  of  several  Tesla.  These  fields  are  achievable 
either  with  pulsed  magnets  or  superconducting  magnets.  A superconducting  magnet 
is  out  of  the  question  for  explosive  MHD  generators  because  it  requires  a cryo- 
genic system,  which  obviously  eliminates  it  from  applications  requiring  a long 
shelf  life  and  little  or  no  preparation.  This  leaves  only  the  pulsed  magnet. 


Drovided  the  magnet  power  can  be  drawn  from  the  MHD  generator  output.  Self- 


excitation would  be  initiated  by  means  of  a permanent  magnet. 


To  apply  the  self -excitation  concept  to  explosive  MHD  power,  one  more 
barrier  problem  must  be  solved.  Since  an  air  core  magnet  is  an  inductor,  MHD 
power  must  be  generated  for  a sufficient  length  of  time  to  allow  buildup  of  the 
magnetic  field.  For  net  energy  output  in  the  kilojoule  range  from  explosive  MHD 
devices,  the  rise  of  the  field  is  a few  milliseconds,  10  to  100  times  longer  than 
the  test  times  reported  in  previous  explosive  MHD  experiments.  It  has  been 
suggested  that  the  explosively  generated  plasma  be  collected  in  a reservoir  for 
the  necessary  period  of  time  needed  for  self-excitation.  Other  suggestions  have 
included  controlled  collapse  of  a metal  tube  to  generate  long  plasma  test  times 
and  the  use  of  solid  phase  reactions  to  heat  the  test  gas.  Whichever  method  is 
selected,  to  show  the  feasibility  of  explosive,  self-excited,  MHD  power  one  must 
demonstrate  two  items:  1)  The  plasma,  once  produced,  will  remain  in  the  necessary 

state  for  sufficient  time  to  achieve  self-excitation,  2)  Self-excited  MHD  power 
is  feasible  with  the  conditions  produced  in  explosively  generated  plasmas. 

Development  of  the  MHD  generator  is  more  expedient  using  shock-tubes  as 
several  experiments  can  be  done  each  day  and  longer  test-times  are  attainable. 

The  test-times  are  1-14  milliseconds  in  the  two  shock-tunnel  facilities  being 
utilized.  These  facilities  are  the  ONR/GE  MHD  Generator  Shock  Tunnel  located 
at  Valley  Forge,  Pa.,  and  the  other  the  NASA  Ames  Research  Center,  Electric  Arc 
Shock  Tunnel  (EAST)  facility  at  Moffat  Field,  California.  The  GE  tunnel  produces 
plasmas  in  the  10,000°K,  several  atmosphere  range  and,  has  a solenoidal,  air  core 
electromagnet  around  the  channel.  This  channel  is  being  used  for  parametric  studies 

of  MHD  generator  operation  in  the  high  temperature,  low  Hall  parameter  equilibrium 

o 

regime.  The  EAST  facility  produces  plasmas  in  the  15  -*  18,000  K,  100  atmosphere 
range.  The  pulsed,  self-excited,  MHD  generator  (PSEMG)  channel  will  be  developed 
to  match  the  EAST  plasma  conditions.  The  EAST  facility  is  appropriate  for  testing 
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the  PSMEG  concept  under  conditions  of  self-excitation.  Based  on  results  of 
these  tests  plus  future  explosive  heating  development  a lightweight  PSMEG 
package  could  be  designed  for  use  with  an  explosive  source. 

Four  channels  have  been  utilized  in  the  experimental  program  thus  far. 
They  are: 

(1)  a supersonic  channel  (previously  used  for  MHD  laser  work) 

(2)  a subsonic  channel  with  a solinoidal  magnet  - GE  I. 

(3)  A high  pressure,  high  temperature  test  channel  - NASA  I. 

(4)  a channel  similar  to  (3)  but  with  a permanent  magnet  to  look  at 
self-excitation  - NASA  II. 

Channel  (1)  has  been  fully  described  in  Reference  1 and  will  be  dwelt 
briefly  upon  in  the  following  section.  A fu! 1 description  of  the  other  three 
channels  will  be  given  in  section  3 along  with  a description  of  the  experimental 
results . 
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2.  THEORETICAL  DISCUSSION  AND  DESIGN  CONSIDERATIONS 
2.1  In t rod act  ion 

In  the  previous  annual  report^the  feasibility  of  the  pulsed  self-excited 
MED  Generator  (PSEMG)  concept  was  investigated  and  it  was  concluded  that  enthalpy 
extraction  of  the  order  of  10%  was  possible  but  that  detailed  consideration  must 
be  given  to  the  transient  interaction  of  generator,  permanent  magnet,  magnet 
coil  and  load  to  achieve  a successful  design.  A key  question  for  the  experiment 
is  electrode  drop,  especially  at  initiation  when  only  the  permanent  magnet  B 
field  is  present. 

Self-excitation  requires  an  initial  field  supplied  by  permanent  magnets 
(PM).  An  electromagnet  coil  (EM)  in  series  with  the  external  load  must  then 
build  up  the  field  to  its  required  value.  Calculation  of  rise  time  and  equivalent 
inductance  must  take  into  account  the  changes  in  flux  contributed  by  both  the 
permanent  magnet,  up  to  its  saturation  level,  and  the  coil.  Hysteresis  loss 
is  not  significant  since  the  PM  is  never  demagnetized,  but  eddy  currents  must 
be  blocked  in  the  actual  design  by  fabricating  the  PM  from  thin  strips  in  the 
manner  of  a transformer  core. 

The  PSEMG  depends  on  the  magnetic  field  supplied  by  a pair  of  permanet 
magnets  to  induce  an  initial  voltage  and  current.  This  current,  passed  through 
an  electromagnet  coil,  augments  the  field  and  makes  possible  a rapid  buildup 
of  voltage  and  current.  To  analyze  the  complete  circuit  it  was  necessary  to 
obtain  realistic  but  mathematically  tractable  models  of  each  component  and 
their  interactions. 
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2_._2 MUD  Generator  Model 

Flow  velocities  are  of  the  order  of  2000  m/sec  and  the  MHD  generator  is 
less  than  20  cm  long,  so  time  of  flow  through  the  generator  is  of  the  order 
of  0.1  milliseconds,  which  is  5%  of  the  nominal  2 ms  test  time.  Thus,  with 
the  possible  exception  of  time  required  for  transition  from  supersonic  to  sub- 
sonic flow  (where  flow  stability  may  be  a consideration),  steady-state  MHD 
generator  operation  is  a reasonable  assumption. 

Under  the  very  high  temperature,  high  pressure,  flow  conditions  appropriate 
to  the  PSEMG,  the  analysis  of  a single  electrode  pair,  Faraday  MHD  generator, 
yields  a set  of  voltage-current  characteristics  for  fixed  values  of  magnetic 

field  distribution.  A quasi-one-dimensiona.1  analysis  of  MHD  flow  of  a semi- 

(1) 

perfect  gas  with  friction  and  heat  transfer  was  used  , with  tabulated  data  for 
electrical  conductivity,  specific  heat  and  speed  of  sound.  Assumption  of  a 
single  electrode  pair  imposes  a boundary  condition  of  constant  voltage, 
consequently  loading  parameter  is  a function  of  axial  location.  Magnetic 
field  was  characterized  by  its  value  at  the  generator  entrance  and  assumed  to 
vary  inversely  as  channel  width  - and  also  inversely  as  channel  volume,  since 
channel  height  was  held  constant. 
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Voltage-current  characteristics  at  design  stagnation  conditions  ere  shown  in 
Figure  2-1  for  the  first  MHD  Channel  (designated  Ames  II)  to  be  tested  at  the  NASA 
Ames  EAST  facility.  They  can  be  idealized  as  linear,  with  a di'scontinuity  and  change 
in  slope  at  the  transition  from  supersonic  to  subsonic  flow.  Complete  characteri- 
zation of  the  .MHD  generator  thus  requires  only  a numerical  tabulation  of  five 
quantities  as  a function  of  magnetic  field: 

(1)  Short  circuit  current 

(2)  Transition  current 

(3)  Open  circuit  voltage 

(4)  Voltage  with  supersonic  flow  and  transition  current 

(5)  Voltage  with  subsonic  flow  and  transition  current 

The  circuit  equivalent  is  simply  an  EMF  in  series  with  a resistance,  where 
both  are  B-field  dependent  and  the  resistance  has  different  values  for  supersonic 
and  subsonic  flow.  As  can  be  seen  from  Figure  2-1  , these  five  quantities  are 
smopth  functions  of  magnetic  field  and,  with  the  exception  of  short  circuit  current, 
are  nearly  straight  lines.  Short  circuit  current  reaches  a peak  and  then  decreases 
for  sufficiently  strong  magnetic  field,  which  is  a reasonable  result  of  very  strong 
magnetic  interaction. 

2 . 3 Electromagnet  Coil  Model 

2 3 

A crescent  coil  model  was  used  for  magnet  coil  design  calculations.  The 
crescent  coil  model  assumes  that  the  coils  are  distributed  as  shown  in  the  sketch. 
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The  size  of  the  coils  is  governed  by  a single  parameter,  the  ratio  C/D^,  and  resulrs 


are  nearly  identical  ti  a rectangular  coil  model  in  which  optimum  combinations  of 
two  geometry  parameters  must  be  determined.  The 
transient  field  resulting  from  a suddenly  applied 
constant  voltage  can  be  computed.  The  effects  of 
end  turns  and  of  a change  in  coil  resistance  due  to 
2 

I R heating  are  included  although  the  latter  is 
negligible  on  the  millisecond  time  scale  of  PSEMG  operation. 

From  a PSEMG  circuit  standpoint,  the  performance  of  the  electromagnet  coil 
is  completely  characterized  by  resistance  and  inductance,  plus  magnetic  flux  density 
as  a function  of  current  to  the  coil.  For  a given  magnet  bore  (MHD  channel  plus 
clearances),  the  ratio  of  ampere-turns  to  magnetic  flux  density,  the  reciprocal  of 
the  ratio  of  volts  per  turn  to  flux  density  and  the  ratio  of  inductance  to  resistance 
(i.e.  the  time  constant  of  the  coil)  are  all  linear  or  nearly  linear  functions  of 
C/D^.  Given  these  three  ratios,  specification  of  C/D^  and  the  total  number  of  turns 
is  sufficient  to  completely  determine  the  electromagnetic  characteristics  of  the 
coil.  These  quantities  are  shown  in  Figure  2-2  for  the  Ames  II  design. 

2 . 4 Magnetic  Field  Model 

Magnetic  flux  density  inside  the  MHD  generator  channel  is  the  sum  of  the  flux 
density  due  to  the  permanent  magnets  and  that  due  to  the  coil.  A current,  i,  in  the 
coil  produces  an  applied  magnetic  intensity,  H , proportional  to  i and  a corresponding 
flux  density,  B=  U0H^.  The  applied  H also  changes  the  value  of  inside  the  permanent 
magnet  material  and,  until  it  saturates,  results  in  an  added  contribution  from  the 
permanent  magnet  as  illustrated  in  Figure  2-3. 

Figure  2-3(a)  shows  conditions  for  no  current  in  the  coil,  H=0.  The  state  of 


the  permanent  magnet  is  given  by  the  intersection  of  the  load  line  which  is  a function 
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of  magnetic  circuit  geometry,  with  the  demagnetization  curve,  which  is  a 
property  of  the  magnetic  material,  see  Appendix  A.  Fringing  reduces  the  flux 
density  in  the  air  gap  to  the  level  indicated.  In  Figure  2-3 (b)  some  current  flows 
in  the  coil,  resulting  in  a non-zero  value  of  H and  a corresponding  flux  density, 
Bc,  due  to  the  coil  of  U0H0.  The  load  line  for  the  permanent  magnet  circuit  now 
originates  at  the  applied  value  of  H and  its  intersection  with  the  demagnetization 
curve  curve  shifts  up  and  to  the  right,  resulting  in  a larger  contribution,  PM, 
from  the  permanent  magnet  in  addition  to  Bc . Once  Che  magnetic  material  reaches 
satur  .tion,  there  is  no  further  increase  in  B as  shown  in  Figure  2-3(c).  Total 
flux  density  is  the  sum  of  3pj.j  and  Bq  as  shown  in  Figure  2-3  (d).  In  Figure  2-3 
the  demagnetization  curve  has  been  idealized  as  two  straight  lines  intersecting 
at  the  point  of  saturation.  For  samarium  cobalt  this  approximation  is  valid  to 
within  a few  percent. 


2 . 5 Circuit  Model  for  Combined  Coil  and  Permanent  Magnet 


To  apply  the  above  behavior  to  an  electro- 
magnetic circuit  model  of  the  PSEMG, 
consider  the  circuit  shown  in  the  sketch. 
Summing  voltages  around  the  loop 


V = iR.  + l -!i-+  v d0PM 
C c dt  dt 


(1) 


Where  V is  net  output  voltage  (generator  voltage  minus  electrode  drop),  Lc  is  the 

inductance  of  the  air  core  coil,  _d  PM.  is  the  rate  of  change  of  flux  produced 

at 

by  the  permanent  magnet  due  to  the  field  produced  by  the  coil  and  linking  the  coil 
and  N is  the  number  of  turns.  In  samarium-cobalt  magnet  material,  magnetic  flux 
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density,  B,  varies  linearly  with  magnetic  intensity,  H,  up  to 

flux  0pm  is  proportional  to  B and  since  H is  proportional  to 

follows  that  ps  pr0port:ional  to  —^2—  and  the  constant 

dt  dt 

is,  in  effect,  an  extra  inductance  which  is  in  the  circuit  on 
occurs.  Equation  1 can  thus  be  rewritten 


saturation.  Since 
current,  i,  it 
of  proportionality 
ly  until  saturation 


V = iRc  + (Lc+M)-7r 


(2) 


Li  = Lpj^  (PM  unsaturated) 

Lj^  = 0 (PM  saturated) 


For  further  details  see  Appendix  B. 

An  additional  back  EMF  resulting 
can  be  avoided  by  laminating  the  magnet 
estimating  eddy  currents  indicated  that 
sufficiently  fine  to  avoid  the  problem. 


from  eddy  currents  induced  in  the  magnet 
material.  Two  difference  models  for 
a lamination  thickness  of  1/16"  was 


2 . 6 PS  EMC.  Circuit 

Figure  2-4  shows  schematically  the  basic  components  of  the  PSEMC,  circuit. 

The  MHD  generator  is  represented  as  a resistance  and  an  EMF,  both  dependent  on  the 

magnetic  t.old  produced  by  the  coil.  A separate  external  EMF  simulates  electrode 

voltage  drop.  Jenerator  output  passes  through  a coil  wtiich  has  a resistance,  R, 

2 

and  an  inductance,  l,f..  Coil  resistance  is  assumed  constant  because  i R heating  is 
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negligible  and  skin  effect  can  be  minimized  by  using  multi  strand-woven  wire  (Litz 
wire).  Inductance,  Lq,  includes  the  effect  of  the  initially  unsaturated  permanet 
magnet.  The  load  is  presumed  to  be  some  combination  of  resistance  and  inductance 
because  of  the  inherently  high  current,  low  voltage  nature  of  the  PSEMG. 

A diode  is  shown  in  parallel  with  the  MHD  generator.  When  the  plasma  flow  is 
exhausted  the  generator  will,  in  principle,  become  an  open  circuit  at  a time  of 
maximum  current  in  the  circuit.  The  diode  allows  the  current  to  flow  after 
generator  cutoff  so  that  inductive  energy  stored  in  the  coil  and  load  is  delivered, 
in  part,  to  the  load  resistance.  For  initial  experiments,  this  function  could  be 
served  by  a protective  spark  gap  or  the  driver  gas  may  be  sufficiently  conductive. 

All  circuit  elements  and  couplings  having  been  reduced  to  linear  or  piece-wise 
linear  models  for  analysis  of  transient  behavior,  they  were  then  incorporated  into 
a time-share  computer  code,  COUP.  For  a given  MHD  generator,  transient  circuit 
performance  was  explored  as  a function  of: 
o Load  Resistance 
o Load  Inductance 
o Electrode  drop 
o Number  of  Turns  in  the  Coil 
o C/D^  (i.e.  Coil  size) 

Initial  voltage  is  determined  by  the  open  circuit  MHD  generator  in  the  field 
of  the  permanent  magnet.  Current  at  any  instant  is  found  by  a time  integration  of 
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the  complete  circuit  equation.  This  in  turn  determines  the  instantaneous  magnet 
system  and  generator  performance.  Total  energy  delivered  to  the  coil  and  to  an 

external  load,  if  any,  is  also  found  by  integration. 

An  upper  bound  calculation  on  energy  delivered  by  the  generator  with  no 
load  except  the  magnet  coil  itself,  Figure  2-5  , illustrates  the  importance 

of  optimization  of  the  magnet  coil  configuration.  For  a given  number  of  turns, 
total  energy  delivered  by  the  generator  is  an  optimum  for  some  value  of  the 
combined  resistive-inductive  load  representing  the  coil  and  4 or  6 turns  best 
matches  the  time  varying  generator  output  voltage.  Only  even  numbers  of  turns 
are  shown  to  maintain  symmetry  with  half  the  coil  above  the  channel  and  half 
below.  At  least  4 turns  are  desirable  to  minimize  the  perturbations  resulting 
from  connecting  leads  to  electrodes  and  to  any  external  load. 

Matching  of  external  load  to  the  generator  and  coil  is  also  important  as 
illustrated  in  Table  2- 1 for  a purely  resistive 'load . For  a 1 milliohm  load, 

4 turns  is  clearly  optimum  and  an  increase  or  decrease  in  load  resistance  de- 
creases the  energy  delivered.  The  locus  of  generator  operation  for  the  1 milli- 
ohm, 4 turn  case  is  plotted  on  Figure  2-  L Note  that  steady-state  is  not  reached 
inr  the  2 ms  flow  time.  While  not  shown  on  Figure  2-1,  energy  stored  in  the  field 
of  the  magnet  is  delivered  to  the  load  for  about  2 ms  following  cutoff. 
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Several  channels  have  been  designed  and  built,  with  applications  discussed  in 
Section  3 of  this  report.  The  latest  channel,  designated  AMES  II  was  the  first 
attempt  at  a self-excited  package  design.  Figure  2-6  .The  channel  has  a throat, 
sized  on  the  basis  of  predicted  flow  time  for  the  NASAnEAST"faci lity  of  2 ms, and 
a 2 to  1 expansion  ratio. 

A single  pair  of  tungsten  coated  copper  electrodes  have  external  leads  to 
permit  series  connection  with  the  magnet  coil  or  separate  excitation  of  the  coil. 

The  permanent  magnets  were  fabricated  by  Hitachi  Magnetics,  Inc.  from  a 
samarium-cobalt  alloy  with  a demagnetization  curve  shown  in  Appendix  A. 

The  magnets  were  fabricated  from  thin  laminations  cemented  together 
with  an  epoxy  resin.  The  laminations  are  sized  to  limit  eddy  current  loss  but 
were  also  necessary  from  a manufacturing  standpoint  since  the  magnets  are  much 
too  large  to  make  as  a single  unit.  A laminated  soft  steel  platform  was  necessary 
for  structural  support. 

The  reluctance  of  the  return  path  is  about  25  times  that  of  the  channel  gap. 
However,  no  return  path  was  provided  because: 

(a)  External  iron  can  "short  circuit"  the  flux  in  the  air  gap. 

(b)  Extra  iron  increases  inductance. 

(c)  There  is  considerable  structural  conflict 

The  entire  unit  was  filament  wound  with  epoxy  impregnated  fibreglas  to 
provide  a structure  capable  of  repeated  shock  tunnel  testing  with  an  ample  factor 
of  safety. 

The  models  discussed  herein  are,  of  course,  idealized  in  many  ways.  Experi- 
mental results  may  indicate  that  modifications  to  some  of  the  assumed  inputs  or  to 
the  models  themselves  may  be  necessary  to  provide  an  accurate  prediction  of  PSEMG 
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performance  and  this  may,  in  turn,  indicate  design  improvements  needed  for  future 
PSF.MG  channels.  For  example,  preliminary  data  using  the  Ames  II  channel  have 
indicated  that  stagnation  pressure  achievable  in  the  EAST  facility  is  less  than 
that  assumed  in  the  Ames  II  design,  as  discussed  in  Section  3. 

The  analysis  has  already  shown  that  a substantial  fraction  of  the  generator 
output  is  devoted  to  self-excitation  although  some  of  the  energy  can  be  recovered 
in  a resistive  load  subsequent  to  generator  cutoff.  If  the  primary  objective  were 
the  creation  of  a magnetic  field,  the  PSEMG  would  be  most  effective  if  the  self- 
excitation field  itself  were  applied  to  the  end  use. 

In  the  developmental  stage,  emphasis  has  been  on  a materials  and  structural 
design  best  suited  to  repetitive  testing  with  the  shock  tunnel.  Structural  integrity 
has  been  the  paramount  consideration  with  no  attempt  to  minimize  size  and  weight  of 
the  overall  package.  Once  operating  conditions  are  well  established,  a considerable 
size,  weight  and  cost  reduction  is  possible. 
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3.  EXPERIMENTS  AND  DESCRIPTION  OK  HARDWARE 


3 . 1 ONR/GE  Shock  Tunnel 

The  ONR/GE  facility,  Figure  3-1  , is  a reflected  shock  tunnel  with  a 6 inch 

diameter  driver  and  a 12  inch  diameter  driven  tube.  It  was  designed  primarily 
to  deliver  plasma  for  testing  up  to  10  ms  at  the  relatively  modest  temperature 
and  pressure  (typically  3500°K,  5 atm)  appropriate  for  a non-equilibrium, noble 
gas  MHD  generator.  The  driver  can  be  pressurized  to  2500  psi  and  is  equipped 
with  external  electric  heaters  to  raise  the  driver  gas  temperature  to  approxi- 
mately 600°K.  These  are  the  upper  bounds  on  driver  energy.  Performance  then 
depends  on  the  gas  in  the  driven  tube  and  the  driven  tube  pressure  p^. 

Predicted  performance  with  argon  is  showTi  in  Figure  3-2  . Figure  3-2 
shows  that  stagnation  temperature  drops  rapidly  above  p^=5  atm  so  an  upper  bound 
reference  condition  was  taken  as  12,000°K  and  5 atm.  Figure  3-2  also  shows 
measured  shock  Mach  numbers  obtained  during  the  preliminary  testing  discussed 

later.  A test  time  of  approximately  1.5  ms  was  obtained  at  the  Mg  = 7.6, 

2 

p^  = 5 torr  condition  with  a channel  having  a 50  cm  throat. 

A capacitor  powered  electromagnet  capable  of  delivering  a uniform  field 
of  more  than  3 Tesla  for  more  than  10  milliseconds  is  available, as  shown  in 
Figure  3-j.  With  this  electromagnet,  any  desired  initial  field  can  be  applied 
without  the  encumbrance  (and  expense)  of  a permanent  magnet  installation.  The 
facility  is  therefore  well  suited  to  exploratory  testing  and  diagnostic 
measurements . 

Initial  experiments  were  done  with  an  available  MHD  channel  to  get  the  test- 
time stagnation  temperature  (T5)  and  stagnation  pressure  (P5)  capability  of  the 
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shock-tunnel . At  a hj  of  5 mm  Hg  and  of  1200  psia  the  test-time,  T5  and  P5 
were  1.5  ms,  10,600°K  and  2 atmospheres  respectively. 


The  MHD  experiments  done  with  this  channel  have  been  fully  reported  in  Reference 
1 and  only  a brief  summary  is  given  here.  The  Hall  parameter  was  less  than  unity 
and  the  channel  was  operated  with  all  electrode  leads  connected  externally  to  give, 
in  effect,  a single  electrode  pair.  The  power  output  dropped  40%  as  compared  to 
segmented  electrode  operation,  which  proved  that  the  experimental  conditions  were 
already  in  a regime  where  segmentation  of  the  electrodes  is  not  critical.  Channel 
construction  is  considerably  simplified  with  the  corollary  benefits  of  simplified 
measurements  and  analysis.  In  the  Ames  II  experiments  the  pressure  was  nearly  two 
orders  of  magnitude  higher  so  the  Hall  effect  was  negligible. 

The  conclusions  from  this  series  of  experiments  were  that  the  shock  tunnel 
was  capable  of  producing  a plasma  with  a conductivity  large  enough  to  do  meaningful 
high  magnetic  Reynolds  numbers  experiments  with  sufficient  test  time  (•>■'1.5  ms). 

The  next  channel  to  be  described  (GE  I)  was  specifically  designed  for  high  magnetic 
Reynolds  number  operation  and  its  geometry  tailored  to  conditions  achievable  in 
the  shock-tube. 

3.2  GE  I MHD  Channel 

The  GE  I channel  was  designed  to  stu^y  MHD  operation  at  high  temperatures. 
Drawings  of  the  assembly  are  shown  in  Figures  3-3  and  3-4  and  a photograph 
of  the  channel  in  place  on  the  shock-tube  in  Figure  3-5.  The  channel  has  a 
throat  size  of  5 cm  x 8 cm  and  has  a divergence  such  as  to  keep  the  flow  sub- 
sonic when  operating  at  designed  power  output.  The  geometry  of  the  channel  was 
described  fully  in  Reference  1.  There  are  five  pairs  of  copper  electrodes  which 
' 
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can  have  individual  loads  or  can  be  externally  connected  to  form  a single  elec- 
trode pair.  The  five  anodes  have  a .005"  sprayed  tungsten  surface  to  reduce 
oxidation  and  erosion. 

The  plasma  conditions  achievable  in  the  GE/ONR  shock-tunnel  are  not 

sufficient  to  give  full  excitation.  A split-solenoid  air-core  magnet,  powered 

by  a 2 farad  capacitor  bank,  is  therefore  used  to  simulate  a permanent  magnet 

field.  This  gives  experimental  flexibility  to  pick  the  magnetic  field  and  also 

» 

freedom  to  explore  the  channel  operation  as  a pure  MHD  generator.  A coil  for 
self-excitation  experiments  is  also  placed  around  the  channel.  It  has  3 turns 
above  and  3 turns  below  the  channel. 

The  operation  of  this  channel  in  a standard  Faraday  MHD  generator  con- 
figuration has  been  explored.  Some  shorting  of  the  current  downstream  at  the 
diaphragm  holder  was  experienced.  Replacement  of  aluminum  flanges  with  fiber- 
glass reinforced  epoxy  and  teflon  spraying  of  other  metal  surfaces  minimized  the 
problem.  Current  at  four  of  the  electrode  pairs  are  shown  on  Figure  3-6  and 
the  voltage-current  characteristics  obtained  are  compared  to  theoretical  calcu- 
lations on  Figure  3-7  . The  theoretical  curves  do  not  include  any  electrode 
or  end  losses.  The  power  outputs  obtained  are  thus  very  good  although  the 
magnitudes  of  the  apparent  voltage  losses  are  higher  than  past  experience  with 
General  Electric's  low  temperature  Faraday  generators  would  have  led  one  to 
expect.  However,  this  generator's  interaction  length  is  only  15  cm  versus  the 
120  cm  of  the  GE  ST-40W  channel  and  thus  end  effects  undoubtedly  play  a much 
greater  role.  Table  3-1  presents  results  of  generator  operation  in  three 
different  configurations  - (a)  simple  MHD  generator  with  5 electrode  pairs 
(b)  simple  MHD  generator  with  a single  electrode  pair  and  (c)  with  the  generator 
output  applied  to  the  augmenting  coil.  Table  3-1  shows  that  the  enthalpy 
extraction  (electrical  power  out/heat  in)  varies  from  0.9%  at  0.5  Tesla  to  11.5% 
at  1.3  Tesla.  The  highest  enthalpy  extraction  obtained  from  the  ST-40W  channel 


was  19%  at  the  much  higher  field  of  2.8  Tesla,  so  in  terms  of  enthalpy  extraction. 


this  generator  is  performing  in  excellent  fashion.  The  power  output  went  down 
during  single  electrode  pair  operation  by  40-50%  compared  to  5 electrode  pair 
operation  (this  is  in  agreement  with  the  results  discussed  previously). 

With  the  augmenting  coil  connected  to  the  electrodes  more  power  (using  the 
same  plasma  conditions)  was  obtained  than  in  the  simple  Faraday  mode  (Table  3-1  ). 

This  was  a combination  of  (a)  better  matching  of  the  generator  impedance  to  the 
external  load,  and  (b)  the  generator  current  applied  to  the  coil  gave  an  effective 
7-»8%  increment  in  B.  A resistor  was  in  series  with  the  coil  for  the  results 
presented  in  Table  II-l  to  give  better  load  matching.  Typical  (Run  36)  voltage 
and  current  oscillograms  are  shown  in  Figure  3-8.  The  average  power  in  35  K watts 
and  the  output  energy  in  56  joules.  This  particular  case  is  for  the  lowest  B field 
used  (0.5  Tesla)  and  as  can  be  seen  from  Table  3-1  the  power  (and  energy)  output 
goes  up  nearly  exactly  as  B squared.  The  emphasis  has  been  on  low  B field  because 
a self-contained  PSEMG  would  be  restricted  to  the  magnetic  fields  obtainable  from 
permanent  magnets.  However  the  power  output  at  higher  B fields  is  of  interest  as 
it  shows  what  power  levels  would  be  reached  if  greater  augmentation  was  achievable. 

So  far  the  power  output  goes  up  as  B sqauared  with  no  levelling  off,  I.e.  end  shorting 
or  any  other  loss  mechanism  present  is  not  becoming  more  important  as  B increases. 

The  conclusions,  thus  far,  from  these  experiments  are  thac  the  MHD  generator 
self-excitation  concept  is  sound  but,  of  course,  plasma  energy  densities  several 
orders  of  magnitude  larger  and  plasma  temperatures  twice  as  large  will  be  required 
to  give  energy  outputs  in  the  kilojoule  range  in  a small  package. 

3.3  NASA  Ames  "EAST"  Facility 

Development  of  the  PSEMG  package  requires  a more  powerful  test  capability. 

The  Electric  Arc  Shock  Tunnel  (EAST),  which  was  built  primarily  for  use  as  a 
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hypersonic  wind  tunnel,  is  admirably  suited.  This  unique  facility,  located  at 
the  NASA  Ames  Research  Center,  Moffett  Field,  California,  is  available  through 
an  operating  agreement  betw-en  DoD  and  NASA.  Testing  at  temperatures  even  higher 
than  the  nominal  17,000°K  will  be  possible. 

Instrumentation  and  experimental  flexibility  are  necessarily  limited  by 
the  requirements  for  tightly  coupling  generator,  coil  and  permanent  magnet  in 
a container  designed  to  withstand  400  atmospheres  safely.  Stagnation  conditions 
and  external  load  are  variable  but  diagnostics  are  limited  to  exit  stagnation 
pressure  plus  electrical  performance. 

3 . 4 Ames  I Channel  Configuration  and  Experiments 

To  check  the  operating  parameters  of  the  NASA  EAST  facility  and  to  check 
channel  construction  techniques  and  materials, a test  channel  (designated  Ames  I) 
was  fabricated.  Figure  3-8  is  a schematic  drawing  of  this  channel.  The  walls 
were  of  boron  nitride  and  there  were  two  copper  electrodes.  The  hannel  diverged 
along  the  electrode  walls  and  was  of  constant  height  (5/b")  between  the  insulator 
walls.  Three  pressure  transducers  were  located  in  the  insulator  walls  to  check 
the  aerodynamic  performance.  The  channel  was  strengthened  to  withstand  the 
expected  high  pressures  by  being  wound  to  an  O.D.  of  8 inches  with  epoxy  im- 
pregnated continuous  wound  "Fiberglas".  Three  different  front  face  plates 
exposed  to  the  high  temperature  were  used  to  check  construction  materials  for 
future  channels.  They  were,  G-10  melamine,  stainless  steel  and  copper.  They 
each  stood  up  to  several  tests,  but  the  copper  face  plate  was  used  for  most  of 
the  tests  . Figure  3-9  is  a photograph  of  the  channel  taken  between  tests  as 
the  channel  is  removed  after  each  run  for  replacement  of  a mylar  diaphragm 
located  at  the  rear  of  the  channel. 

The  EAST  facility  has  three  drivers:  CONICAL,  54"  cylindrical  and  30" 

cylindrical.  It  is  driven  by  a capacitor  bank  which  can  be  connected  in  various 

18 


ways,  two  of  which  are  designated  "20  KV  MODE"  and  40"  KV  MODE".  Each  combination 
of  driver  and  capacitor  bank  configuration  has  its  own  desirable  and  undesirable 
features.  For  instance,  the  30"  cylindrical  driver  will  have  a higher  energy 
density  from  a given  capacitor  bank  configuration  and  operating  voltage  than  the 
54"  driver.  This  results  in  higher  driver  pressure  and  hence  shock  velocity,  but 

test-time  usually  is  shorter  with  a shorter  driver.  An  important  part  of  the  tests 

was  the  determination  of  the  optimum  driver  and  capacitor  bank  configuration  to 
give  large  P^,  and  test-time.  Table  -2  is  a summary  of  the  tests  with 
representative  P^'s,  T^'s,  and  test-times  and  P5  is  shown  plotted  against  Mach 
number  in  Figure  3-10.  The  initial  shock-tube  pressure  (Pj)  was  varied  from  10 
to  50  mm  Hg  and  as  Pj  was  decreased  the  Mach  number  increased  such  as  to  maintain 
P^  fairly  constant.  Operation  with  P^  in  the  20  30  mm  Hg  range  gave  best 

"tailoring"  (i.e.  a flat  P5  pressure  versus  time  profile^.  The  54"  cylindrical 
driver  and  40  KV  bank  operating  mode  gave  conditions  nearest  to  those  hoped  for. 
The  P^'s  were  lower  ('V100  atmos.)  than  hoped  for  400  atmos.)  and  this  lowers 
the  energy  density  of  the  plasma.  The  30"  cylindrical  driver  did  give  higher 
P5’s  but  the  test-times  were  smaller  (1.0  -»1.2  ms)  than  with  the  54"  driver 

(^1.5  ms).  The  P^  shown  in  Table  3-2  for  Run  31  was  a peak  pressure  which 

fell  off  quickly  during  the  test-time.  Representative  oscillograms  obtained 
using  the  54"  driver  are  shown  in  Figure  3-11  . Figure  3-11  (a)  shows  a pressure 
transducer  signal  located  1"  from  the  front  face  of  the  channel.  The  test-time 
was  measured  by  applying  a small  voltage  from  a capacitor-resistor  power  supply 
to  the  channel  electrodes  and  measuring  the  resultant  current;  this  is  shown 
in  Figure  3-11  (b) . An  important  result  of  this  test  was  that  the  plasma  broke 
down  with  less  than  1 volt  applied  to  the  electrodes.  Thus,  the  total  electrode 
voltage  loss  is  very  small  in  the  absence  of  a magnetic  field.  The  test-time 
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was  further  verified  by  a photomultiplier  at  the  same  axial  location  as  the 
P5  measurement.  The  photomultiplier  had  a 7000A  filter  to  look  at  two  argon 
lines  (6965  and  7067A) . A photomultiplier  output  is  shown  in  Figure  3-11  (c) . 

The  conclusion  from  this  series  of  tests  were  that  the  NASA  EAST  facility 
is  capable  of  producing  the  high  conductivity  plasmas -necessary  for  MHD  self- 
excitation experiments. 

t . 5 Ames  1 1 Channel  Configurat  ion  and  Experimen t s 

The  Ames  II  channel  is  similar  to  Ames  1 in  internal  channel  geometry  but 
with  a larger  outlet  to  inlet  area  ratio  to  offset  the  gas  deceleration  as  a re- 
sult of  J x B forces.  The  entrance  area  is  3.5  cm  x 1 cm  and  the  exit  area  is 
7.0  cm  x 1 cm.  The  overall  length  is  16  cm.  The  inside  surfaces  are  boron 
nitride  with  two  copper  electrodes.  The  internal  boron  nitride  and  electrode 
assembly  is  shown  in  Figure  3-12  (TOP).  This  assembly  fits  inside  the  permanent 
magnet  assembly.  Figure  3-12  (Center).  The  permanent  magnet  is  fabricated 
from  samarium-cobalt  and  is  laminated  (with  epoxy  bonding)  to  reduce  eddy 
current  losses.  This  magnet  provides  a field  in  the  active  region  of  the 
channel  uniform  to  within  2%.  A field  of  between  0.2  and  0.3  Tesla  had  been 
hoped  for  but  the  necessary  epoxy  bonding  and  supporting  structure  reduced  the 
gap  field  to  0.14  Tesla.  There  are  two  magnets,  one  above  and  one  below  the 
channel . 

A coil  having  2 turns  wrapped  around  the  top  permanent  magnet  and  2 turns 
around  the  bottom  permanent  magnet  is  connected  to  the  channel  electrodes 
to  provide  the  self-excitation  magnetic  field.  The  coil  was  fabricated  using 
Litz  wire  to  reduce  its  AC  resistance,  Figure  3-12  (bottom).  Holes  are  drilled 
in  the  wrapping  for  connection  of  the  coil  and  electrodes  and  for  instrumentation. 

Characterization  of  this  channel  is  being  done  on  the  EAST  facility  and 
experiments  with  this  channel  have  just  been  initiated.  The  first  experiment 


20 


was  done  with  the  electrodes  open  circuited  and  the  UBD  voltage  measured  - 
Figure  -1-1 3.  The  voltage  was  correct,  showing  that  nothing  unusual  was  happening 
to  the  flow.  In  bootstrap  operation  (electrodes  connected  to  the  augmenting 
coil)  the  power  output  was  very  small  - several  amps  at  several  volts.  The 
channel  was  then  operated  as  a simple  MHD  generator  using  the  permanent  magnet 
to  provide  the  field  and  with  the  augmenting  coil  disconnected.  The  generator 
performed  as  if  its  internal  impedance  were  several  hundred  milliohms  rather 
than  the  3-4  milliohms  expected  from  an  argon  plasma  at  100  atmospheres  and 
17  - 19,000°K.  The  power  output  rose  sharply  for  20  - 30  ys  and  then  dropped 
abruptly.  There  are  several  reasons  which  could  possibly  cause  this  - end 
shorting,  electrode  voltage  losses  at  high  currents,  the  small  initial  magnetic 
field  (.14  Tesla)  and  consequent  small  voltage  to  overcome  losses  and  black-body 
radiation  losses.  However,  the  prime  suspect  is  the  residual  carbon  deposits 
which  are  vaporized  off  insulating  materials  in  the  EAST  driver  by  the  driver 
arc  (300-500  K amps).  As  mentioned  previously  a 1 volt  scource  applied  to  the 
electrodes  was  sufficient  to  break  down  the  plasma,  this  scource  could  only  supply 
several  amps.  A new  low  impedance  source  (40  capacitors  each  of  500  uf  capacity) 
was  constructed.  The  charge  was  sufficient  for  supplying  1,000  amps  for  several 
milliseconds  in  order  to  check  if  the  electrode  losses  increased  at  high  current 
density.  A typical  current  oscillogram  is  shown  in  Figure  3-14  and  it  can  be  seen 
that  the  current  again  starts  rising  quickly  and  then  drops  sharply  after  only 
50  psec.  Again  it  is  suspected  that  losses  are  being  incurred  as  a result  of 
carbon  from  the  shock-tube  driver.  To  check  if  this  is  so,  several  tests  are  going 
to  be  done  as  delineated  in  the  next  section. 
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4.  CURRENT  STATUS 


I 

To  check  if  carbon  is  the  major  loss  factor,  the  EAST  facility  will  be 
operated  pressure  driven.  The  energy  density  will  be  low  but  if  a clean  uniform 
plasma  exists  then  current  (whether  applied  or  generated)  should  be  present  for 
the  full  test-time.  If  this  proves  to  be  the  case,  then  steps  will  be  taken  to 
ensure  cleaner  driver  operation  in  the  arc  mode.  The  present  driver  insulating 
liner  is  RTV  rubber  coated  melamine  and  the  RTV  is  the  primary  source  of  carbon. 
This  liner  will  be  replaced  by  a "Delrin"  liner  already  being  fabricated  - this 
type  of  liner  has  been  used  in  GE  arc-driven  shock-tubes  to  provide  cleaner 
operation 

Another  potential  source  of  difficulty  is  the  low  UBD  voltage  provided  by 
the  permanent  magnet  - approximately  15  volts  open-circuit  are  generated  by  the 
.14  Tesla  permanent  magnet.  To  check  the  generator  operation  at  larger  induced 
voltages  current  will  be  applied  to  the  bottstrap  coil  to  provide  a field  of 
approximately  1.0  Tesla.  A capacitor  bank  to  provide  the  coil  current  is  presently 
being  constructed.  It  will  give  about  30  K Amps  with  a 300  p sec  risetime  and 
will  be  operated  in  a crow-bar  mode. 
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Figure  2-1.  MHD  Generator  Characteristics: 

Single  electrode  pair,  inlet  area  3.5  cm^ , area  ratio  2:1, 
To=17000°K,  Po=400  atm,  132  Kilojoules  thermal  energy  in  2 ms  flow 
time.  Also  shown  is  locus  of  PSEMG  performance  for  1 milliohm 
resistive  load,  0.2  Tesla  PM  field,  20  volt  electrode  dropj  PM 
inductance  included.  Insets  show  key  values  of  steady  state  MED 
generator  voltage  and  current  as  a function  of  magnetic  field. 
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(a)  No  current  in  coil,  H=0 


Saturation 


PM  Demagnetization  Curve 
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(b)  Current  in  coil,  H>0,  PM  not  saturated. 


Figure  2-3.  Magnetic  Flux  density  in  the  MHD  Channel  resulting  from  the 
combination  of  permanent  magnets  and  electromagnet  coil. 
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Figure  2-5 


Upper  bound  calculation  of  energy  delivery  by  PSEMG 
with  magnet  coil  as  the  only  external  load. 

PM  field  = 0.5  T,  electrode  drop  and  PM  inductance 
neglected.  MHD  generator  inlet  area  4.84  cm^,  TQ=17 , 000°K, 
Po=400  atm,  183KJ  total  energy  in  2 ms  flow  time. 


Figure  2 6.  Ames  II  Channel  Design 
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FIGURE  3-7.  VOLTAGE  - CURRENT  CHARACTERISTICS  OF  GE  I MUD  CHANNEL 


STAGNATION  PRESSURE  - PS  I 
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FIGURE  3-10.  MACH  NUMBER  AND  STAGNATION  PRESSURE  OBTAINED  USING  VARIOUS 
EAST  DRIVER  CONFIGURATIONS 


mnImm 


Figure  3-14.  Current  through  plasma  using  a low  impedance 
(high  current)  service 
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Load  Coil  MHD  Energy  To 

Resistance  Turns  Energy  Load 


Milliohras  - Kilojoules  Kilojoules 


1.0  2 1.03  .82 

1.0  4 3.12  1.54 

1.0  6 1.32  .39 

0.5  4 3.81  1.22 

2.0  4 1.83  1.22 


Table  2-1. 


Matching  of  magnet  coil  and  resistive  load; 

C/D  £ = . 08,  20  volt  electrode  drop,  2 ms  flow  time. 
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RUN 

B FIELD 
(TESLA) 

NO.  OF 
ELECTRODES 

rL 

(mii) 

I (AMPS) 

V (VOLTS) 

POWER(KW) 

% ENTHALPY 
EXTRACTION 

30 

0.5 

1 

100 

400 

40 

16 

1.0 

25 

II 

1 

46 

600 

25 

15 

0.93 

28 

II 

1 

22 

800 

18 

14 

0.87 

17 

II 

5 

46 

725 

33 

24 

1.5 

30 

1.0 

1 

100 

1125 

70 

49 

3.1 

25 

II 

1 

46 

600 

52 

58 

3.6 

28 

II 

1 

22 

800 

33 

49 

3.1 

15 

II 

5 

46 

1450 

67 

97 

6.1 

18 

1.35 

5 

46 

2000 

92 

185 

11.6 

TABLE  3-1 (a)  PERFORMANCE  OF  CE  I CHANNEL  OPERATING  AS  A SIMPLE 
MHD  GENERATOR 


RUN 

B FIELD 

LOAD 

ENERGY 

AVERAGE 

— 

(TESLA) 

(COIL  + 

xmfl) 

OUT  + (JOULES) 

POWER  (K1 

33 

0.5 

COIL  + 

0 

34 

22 

39 

II 

COIL  + 

2 

41 

32 

36 

II 

COIL  + 

3 

56 

36.5 

41 

II 

COIL  + 

10 

61 

44.5 

34 

1.0 

COIL  + 

0 

160 

114 

38 

It 

COIL  + 

2 

165 

120 

37 

II 

COIL  + 

3 

214 

153 

42 

II 

COIL  + 

10 

183 

131 

35 

1.35 

COIL  + 

0 

228 

163 

TABLE  3-1 (b)  PERFORMANCE  OF  GE  I CHANNEL  OPERATING  WITH  THE  POWER 
OUTPUT  DELIVERED  TO  THE  AUGMENTING  COIL 
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Table  3-2. 

Plasma  Parameters  Obtained  for  Representative  Runs 
With  the  Ames  I Channel 


APPENDIX  A 


CALCULATION  OF  AIR  CAP  FLUX 
PROVIDED  BY  THE  PERMANENT  MAGNETS 


The  PSEMG  has  permanent  magnets  mounted  above  and  below  the 
MHD  channel,  Figure  A-l,  to  provide  the  initial  magnetic  field  needed 
for  self-excitation.  The  gap  between  the  magnets  makes  a substantial 
demagnetizing  force  inevitable,  even  if  a soft  iron  return  path 
is  provided  between  the  external  poles,  and  under  these  circumstances 
the  samarium-cobalt  type  of  permanent  magnet  is  far  superior  to  other 
magnet  materials.  Calculation  of  magnetic  flux  density  in  the  MHD 

★ 

channel  made  use  of  design  methods  adapted  from  Parker  and  Studders 
and  reference  to  equations  and  figures  marked  (P&S)  refer  to  that 
text.  Equations  have  been  modified  by  a factor  of  when  appropri- 
ate so  that  all  units  are  rationalized 


MKS . 


Two  basic  magnet  configurations 
discussed  by  Parker  and  Studders  are 
presented  here.  They  will  be  used  as 
building  blocks  for  the  complete  mag- 
netic circuit  and  also  serve  to  clarify 
assumptions  and  differences  in  formula- 
tion. A magnetic  circuit,  such  as  that 
shown  in  the  sketch,  is  modelled  in  terms  of  flux,  magnetomotive 
force  and  reluctance  (analagous  to  current,  EMF  and  resistance  in 


Parker,  R.  J.  and  Studders,  R.  J.,  Permanent  Magnets  and  Their 
Application , Wiley,  N.Y.  (1962). 


A-l 


an  electrical  circuit) 


$ = BA  = 


mmf  _ HMLM 

*1  = /£  ' 


Neglecting  fringing  for  the  present,  reluctance 


/e  - 


Equations  1 and  2 can  be  combined  to  form  the  dimensionless  quantity 
called  the  demagnetizing  coefficient,  which  is  a function  of  the 
geometry  of  the  magnetic  circuit  configuration. 


u H..  L 
o M g 

Equation  3 is  a condition  which  is  imposed  on  the  permanent  magnet 
by  circuit  geometry.  Graphically  it  is  represented  on  a B vs.  uqH 
plot  by  a line  through  the  origin  with  slope  — {B/uoHm) , negative 
because  is  negative.  This  line  is  called  the  load  line  and  its 
intersection  with  the  demagnetization  curve  for  the  magnet  material 
determines  the  state  of  the  permanent  magnet  in  a specific  magnetic 
circuit  configuration,  Figure  A-2. 

Now  consider  a bar  magnet  in  air.  Parker  and  Studders 
(P&S  pp  164-165)  model  the  reluctance  of  the  external  return  path  as 

two  spherical  poles  each  radiating  out- 

/ ^ ward  to  infinity.  The  following  deriva- 

/ tion  parallels  Parker  and  Studders  except 

/ r\  ^ y\  \ 

/ / V y \ for  the  use  of  rationalized  MRS  units. 

/ 1 For  one  pole,  of  radius,  r,  equation  (2) 

l for  reluctance  becomes  an  integral. 


/i  - / 

r u 

o 


A 


-.8  -.6  -.4  -.2  0 

UqH,  Tesla 


J I 1 1 1 

800  600  400  200  0 

H,  Kiloampere- turns/meter 

Figure  A-2.  Demagnetization  Curve  for  Samarium-Cobalt 
Magnet  Material  and  Load  Line  for 
Configuration  of  Figure  A-l 
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B,  Tesla 
. 4 

. 2 


A-  4 


to  account  for 


An  "equivalent  r 


terms  of  pole  surface  area 


So  that  reluctance  for  any  reasonable  pole  shape  can  be  written 


Let  the  bar  magnet  be  cylindri 


cal  in  shape  as  shown  in  the  sketch 


Combining  equations  1,  8,  and  10,  the 


demagnetization  coefficient  becomes 


small  correction  for 


For  the  geometry  of  Figure  A-l 
L , = 1.73  cm 

gi 

L.,  = 3.  18  cm 

M 

A =104. 32  cm2 

s = 104.32  + 71.93  = 176.25  cm2 

So  that 


■ H.. 
o M 


= 1.03 


The  value  of  B at  the  operating  point  is  the  flux  density  inside 
the  permanent  magnets.  Some  of  the  flux  will  bypass  the  air  gap 
through  fringe  paths.  To  obtain  an  estimate  of  actual  flux  density 
in  the  air  gap,  consider  the  gap  itself  and  the  fringe  paths  as  two 
reluctances  in  parallel.  Parker  and  Studders  give  an  expression 

(P&S  Eq  (3)  pg  158)  equivalent  to  the 
following  for  reluctance  of  fringe  paths 
originations  from  flat  surfaces 


O'cA’f  ' 5 


£n 


x + i'x^  + xL, 


(15) 


1 + 2 


The  width,  w,  of  the  surface  is  the 
perimeter  of  the  magnets.  Let  the  extent, 
x,  of  the  fringing  be  LM/2  (consistent  with 
the  assumption  that  the  remaining  peri- 
meter surface  is  part  of  the  radiating  pole).  For  the  gap  itself, 


from  equation  2 


L 


For  our  gap  configuration  this  formula  checks  the  more  complete  rela- 
tion given  graphically  by  Parker  and  Studders  (P&S  Fig.  5-28  pg  183). 
Then  with: 

x = .64  + 1.59  = 2.23  cm 
W = 45.2  cm 
L = 1 . 7 3 cm 

g 

A =104. 32  cm2 
The  reluctances  are  : 

(U  H ) c = . 03568  cm-1;  (u  /£  ) = .0166  cm"1 

of  o g 

Assuming  that  the  reduction  in  flux  density  in  the  gas  will  be 
proportional  to  the  ratio  of  gap  flux  to  total  flux,  the  following 
correction  on  B can  be  calculated 


B 

_2  = 
B 

P 


i/uo/?)f  + 17 v?  >g 


. 678 


Thus  the  flux  density  in  the  air  gap  should  be  . 37  x . 68  = . 25  Tesla. 

In  the  initial  calculation  of  demagnetization  coefficient, 
the  reluctance  of  the  fringe  path  was  not  included,  so  that  B/u  H 
should  be  slightly  conservative.  However,  the  fringe  path  calculation 
above  did  not  include  the  additional  gap  area  outside  the  useful  gap 
which  results  from  the  extended  magnet  support  platform.  A rough 
correction  for  this  effect  might  be  to  take  the  ratio  of  magnet  area 
to  platform  area  instead  of  the  above  fringing  correction.  Then 


Bg  _ 104.2 

B 19. 37  x 12. 38 
P 

B = . 37  x . 44  = . 16 

g 
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APPENDIX  B 


MODELS  OF  MAGNETIC  FIELD  AND  INDUCTANCE  DUE 
TO  PRESENCE  OF  PERMANENT  MAGNET 

R. 1 Nomenclature 

Ap^  Area  of  permanent  magnet  pole  faces. 

B^  Magnetic  flux  density  due  to  coil. 

Bp^  Magnetic  flux  density  due  to  permanent  magnets. 

BpMS  Saturation  magnetic  flux  density  in  the  permanent  magnets. 

B Residual  induction.  Flux  density  in  permanent  magnet 

when  Hp^  = 0. 

£p^  EMF  induced  in  PSEMG  circuit  by  changes  in  flux  density  due  to 
the  permanent  magnets. 

H Magnetic  intensity  due  to  coil. 

H Absolute  value  of  intrinsic  coercive  force.  Value  of  H ,,  when 

cPM  _ n PM 

BPM  ' °- 

H Magnetic  intensity  in  permanent  magnet  material  (negative  when 

PM  H = 0.) 

c 

H„,,„  Value  ofHPMwhen  B„.,  reaches  saturation. 

PMS  r‘l  PM 

i Current  in  PSEMG  coil. 

l.p  Inductance  due  to  permanent  magnet. 

1 Number  of  turns  in  magnet  coil. 

s Absolute  value  of  slope  of  load  line. 

s Slope  of  demagnetization  curve, 

m 

t Time. 

u Permeability  of  free  space,  u = 4tt  x 10  ^ henrys/m. 

o o 

^PM  Flux  due  to  permanent  magnets. 

B. 2 Magnetic  Field 

For  samarium-cobalt  the  demagnetization  curve  is  well  represented  by  two 
straight  lines  as  shown  in  Figure  B.l,  so  that 


B-l 


Figure  B.l.  Magnetic  Flux  Density  Due  to 
Permanent  Magnets 


Vt  " "W’pM  + ,JolicPM)  HPM  - HPMS 


R = R 
PM  PMS 


H > LI 

PM  - PMS 


where 


Sm  BRPM^UoHcPM^ 

The  load  line  (see  Appendix  A)  for  the  permanent  magnet  system  is  given  by 
B„ 

where 


PM 


s (u  H - u H ) 
«.  o c o PM 


Sl  BPM7  (’-'oHPM^  1 

The  operating  point  of  the  permanent  magnet  is  the  intersection  of  the 
load  line  and  the  demagnetization  curve.  Eliminating  H between  equations 
1 and  4,  the  PM  flux  density  is 


PM 


s0s  (pi  H + M H ) 
Jc  m o c o cPM 

s + s 
i m 


bpm  bpms 


(hpm  1 W 


(hpm  i hpms) 


(la) 

(lb) 

(2) 

(3) 

(4) 


(5a) 

(5b) 


The  magnetic  field  in  the  MHD  channel  is  the  sum  of  the  field  due  to  the  coil 
and  that  due  to  the  permanent  magnets,  with  the  latter  reduced  by  air  gap 
fringing  (Appendix  A). 

B.  3 Inductance 

Changes  in  magnetic  flux  density  due  to  the  permanent  magnets  will  result 
in  an  EMF  in  the  PSEMG  circuit  given  by 


d0 


£PM  N dt 


PM 


dB 


= NA_ 


PM 


PM  dt 

This  effect  can  be  treated  as  an  additional  inductance  which  vanishes  when  the 
permanent  magnets  are  saturated.  All  the  PM  flux  is  assumed  to  link  the  coil 
so  for  the  purpose  of  calculating  inductance  there  is  no  reduction  for  air  gap 
fringing. 


(6) 


B-3 


J 


dRP.M  = / Vm_\ 

dt  \sc+sm  / 


B.  = 

U H 

l 

O i 

dB 

c 

dt 

(B 

H and  noting  that  i,  H is  a constant, 
c o c PM 


PM  - PMS 


(7) 


But,  for  a given  coil,  B is  proportional  to  the  current  in  the  coil  so  that 


£PM  LPM  -7—  ~ ^ APM 
dt 


(s  s \ , d . 

m | i 1 

s +s  / (i/B)  dt 

V.  m ' r 


(BPM  - bpms^ 


(8) 


The  above  results  havg  been  incorporated  in  code  COUP  for  computing  PSEMG 
performance.  Appendix  C. 
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APPENDIX  C 


COMPUTER  PROGRAM  COUP  AND  SUBROUTINES 


C. I Int  roduct ion 

Code  COUP  and  its  subroutines  were  developed  to  perform  the  time  dependent 
analysis  of  the  PSEMG,  using  models  described  in  Appendix  B and  in  the  main  body 
of  this  report.  Details  of  the  generator  performance  map  are  contained  in  the 

•k 

previous  annual  report.  Time  dependent  behavior  of  voltage,  current,  magnetic 
field  and  output  power  are  computed,  along  with  total  energy  delivered  to  the 
magnet  and  to  an  external  load.  The  integration  uses  a fourth  order  Runge-Kutta 

technique. 

Section  C.2  describes  the  inputs  required.  Section  C.3  contains  sample  runs 
with  explanatory  notes.  Section  C.4  is  a listing  of  the  code  and  its  subroutines 
(except  standard  library  routines  used  for  interpolation  and  integration).  Sec- 
tion C.5  is  a complete  dictionary  including  subroutine  names. 

All  dimensional  quantities  are  in  rationalized  MRS  units  except  for  times, 
which  are  input  in  milliseconds  and  converted  internally. 

C. 2 Input  Data 

First  Level: 

TMAX  Maximum  run  time,  milliseconds. 

TCUT  Time  of  cutoff  of  MHD  power,  milliseconds. 

DTPO  Time  interval  for  intermediate  printout,  millisecond;  . 

Second  Level : 


RL  Load  resistance,  ohms. 

LL  Load  inductance,  henrys. 

VE  Electrode  voltage  drop,  volts. 


* 


Marston,  C.  H.  , Tate,  E. , and  Zauderer,  B.,  "Annual  Report,  MHD  Generator 
Investigations,"  Contract  N00014-73-C-0039,  Project  Code  9000,  1975. 


Third  Level: 


N Total  number  of  turns  in  magnet  coil. 

C C/p^  ^°r  e1uiva^e,nt:  rrescent  coil  magnet. 

LPR  Switch  for  intermediate  printout.  0 = end  results 

only.  1 = print  at  time  intervals  specified  by  UTRO. 

Inputs  are  on  three  levels  to  facilitate  the  search  for  oDtimum  loading  and 
optimum  coil  configuration.  Sett'ng  N = 0 returns  program  for  new  values  of  RL, 
LL  and  VE.  A negative  value  for  RL  or  N returns  program  for  new  time  inputs. 

C . 3 Sample  Runs 

Two  sample  runs  are  shown.  For  the  for  first  the  magnet  coil  is  the  only 
load  on  the  generator.  For  the  second,  an  electrode  drop  of  20  volts  and  a 
1 milliohm  external  load  are  assumed. 
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Ou 

0. 

0. 

\ . )00 

J. 

2 . S2Jh 

04 

1 . / lot 

Oo 

J. 

o . 

.. . uj,  b rt 

o . ^ 

J. 

a/o  (pX  ^e^rta/  /oqo/j 

^ --  o) 

ZZ  »6> 

Zero  brr  nr i*S  i/aii/f  S'  o f /?lt  LE,  VE 


\oT^r/ 


for  j rJ//r 
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« EL,  v c i &,//// jiff,  rrr/Tiii*'  ' o*c/ 

2^  /dr  f/rc-d'^/c  fd/jP 

=4  ’ I I 


: \ , Li'l  , r'l.lL 

/.dddn-04  l.dlbE-uo  o./d4b-J/ 
1',ijC  = 0 . I 2269bd6b-o2 


i ( 

■id ) 

AJ(  vJLib) 

I ( n.lbd) 

4(1) 

rbbi)  ( .< ) 

H 

.-U  ( 1 1 A i 1 

U ) 

J 

. 

j. 1 d9E 

01 

u • 

2 .ul 2L-01 

O . 

u . 

J 

• 2uo 

d . 9 ddE 

ol 

1 . 2o  1 u 

Od 

2. / 7oE-0 I 

2 .474  b 

04 

1 .o4uE 

VJO 

J 

• 4 oo 

d. 1 22b 

ul 

d.  1 v/E 

Od 

d.92uE-0l 

9 . 9 o 1 c 

u4 

1 . u22b 

o4 

J 

. OUO 

o .OO  YE 

ol 

o .02db 

Od 

b .oO  /E-j 1 

2 . / / )t 

U'J 

4 . Odun 

U^t 

0 

• duo 

7 .dd9E 

ul 

y . 42  / E 

Od 

/. o9  /E-u l 

j . d«2  b 

0‘J 

y . u i /E 

1 

• 0 dO 

0 . d<:  on 

ul 

1 . 4ooE 

04 

1 • U 1 Ob  UU 

9 . yybE 

J'j 

1 .oo4E 

OO 

1 

. JdO 

1 . o /dE 

u2 

1 .dd7E 

04 

1 .29  /n  Ou 

1 • 0 o d b 

OO 

d . d /oE 

UO 

1 

.400 

1 .21 VE 

02 

2 . doob 

04 

1 .oOyi:  00 

2 . 4ud  b 

OO 

•j.  bOdE 

UO 

1 

. ">  00 

1 . 4boE 

02 

2 . 229b 

04 

1 . y4vE  OU 

4.474b 

UO 

d . doOE 

UO 

i ) J\  i 

1 . 

/ d je- 

04 

1 

.4  00 

1.4/Ob 

02 

d.  b /oE 

04 

2 . d 1 9E  JO 

4 . do4b 

OO 

1 .2 /9b 

0O 

ljl  J , Ed rLO  , EL0 AU  , Ebb 


J. 2*4 4b  04  I .0d2E  Od  d.4/dE  )2  2.o7bL-Ol 


2.000 

0. 

4.2d4E 

u4 

2 . o44E  00 

0. 

1 .ol uE 

Od 

2.200 

0. 

d. 4d9E 

04 

2.224E  OU 

0. 

1 . 1 49E 

00 

2.400 

0. 

2 . 4 o J E 

04 

1 . y i 9E  00 

0. 

d.292E 

Od 

2.600 

0. 

2 .446b 

04 

1 . oo 1 E Ob 

0. 

o. 9db E 

Od 

2.4  Oo 

0. 

2.  ) 7 d E 

04 

1.441k  00 

). 

4 . d20E 

Ud 

4 .000 

0. 

1 . 1 66E 

04 

1 . 2SbE  Ou 

0. 

d.  1 1 dE 

Ud 

4 . 2 Ou 

0 . 

1 .booE 

04 

1 .uOoE  00 

0. 

2 . 2d  1 E 

oo 

4.400 

0. 

1 .2  lot 

04 

y.oi  /E-o i 

J. 

1 . 02 dE 

Ud 

4 .OOu 

0. 

1 . Jd4b 

04 

d . 47db-0 1 

0. 

1 . 1 /4E 

Uj 

4 .4 uO 

0. 

9.2O0E 

Od 

/.bO lb-01 

0. 

o . ^o4E 

u4 

4.000 

0. 

/ .dloE 

04 

u . o / Ob- 0 1 

0. 

0. 109b 

04 

u xJ, Erb 

1 . d04E 

04 

b . 4 /9E-0 1 

n,0,  LbE 
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C. 4 Program  Listings 


L I SIM  CUU  H ' 

I I / )-_>/  /o  I / .04 


IJ  * JoOHi^EJ  MHL)  ob.4,  ,.iAoi4E)  A.4i)  ._0Ao,  I i.CLUuI  i«J 
2J  : A 1 ,-rt  .4aL  OcKiV 
J J P E AL  I » i Ji« . i- , L,.\ , !.L 
4 ) > , I OC  , A J ) 

OJ  C0MmJ.<  i , Jl  ( j)  , jl  j)  ,L,I( 

0 j I f JUf  i/f  .4  I I Om  . L.»l  . P M 
Ij  j.  , ^ , Vu,  i'CUi  , LL,  PL 

.JJ  . I SC  , \E , SMAA  , LSCls. 

VJ  i Jo , a PH M,  rtC H...  , JMao , C'JioT , B'JHuo,  LH  i 5 

1 Jo  u,  Joti,  H.1L,  dP'-UX,  HMo  i~t 
Mo  ex,  bn 1 CUR 

12o  JIMLhSIJ.4  iE\H(2o) 

I 3o  EQUIVALENCE  ( I ,✓(  1 ) ) , ( EUEh,V(2  ))  , (ELLao,  V(o) ) 
I 4u  EPS® I .E-d 
loo  r F = . 01 

1 >o  14  \i  = o 

I /u  * 

I -i  J * haIJi  UciJ  I”  I r I up 

1 Jo  CAuu  Jb.4 

2 )u  * 

2 1 o * 1.4  HU  T L)a  i A 

22  o I Jo  HP  I IT  s 11 1 4nX , i'CUT , l)  i H(J  (MS)" 

2 Jo  kE/uJii  4AX,  i'CUT,  JiVO 

2 4 o i Maa=TMAX* I . E_3 
2 So  ECU  i'=TCJ  i * I . E—  J 
200  JT  HU=l) i'HO*  I . E-  J 
2/o  I r( i MAX .Li . J. ) SiUH 
2)o  I 10  HP  I IT*  “RL, LL , v E" 

2vo  P E ao  * P _ , LL  , *E 
J )o  IE  (Pc.. Li  . o. ) o J i o I 00 
Jlo  I 2o  HP  1.4  i':  "U.C.LhP" 

J2o  keaU  * .4 , C , L HP 
JJo  lr  (14)  I JJ,  I 10,  I 30 
J4o  1 .3  J CJ.4r  i.40t 
joj  * 

* ■ (IZE  4 /\o.. E a ruP  SPEC  if  itu  H a. 40  C 

3/o  tnu.  AAo 
J J o HP  1.4  T*  "PM , L M , HML  " 

JOo  pP  1 .4 1 J , P.». , LM , H I L 
4 J j H- P.rt+PL 
4 I o L=LM+Lo 
42 o i AuC®( L+HML) /h 
4 Jo  uT*iAUC*fl 

44U  Ir(L)l.jr.l  .E-o)  o 

4-jO  HP  i.4  i' « " 1'aUC®  " , TAuo 


C- 


S 


Ho  J M A. 


1 


460  * 

470  * initialization 

48  O I = 0 . 

4 Oj  b0Ei4=0. 

50u  ELJaU=0. 

510  i =0 . 

520  L0GN=0 
53U  LCUi'=0 
540  TPG=-EPS 

call  RKPd ( uhr i v , temp , r ,dt  , v , u,  w v ) 

5)0  I r (i_PR  .EQ.  U)  00  TG  200 

o/o  PRIui'i!'  T (MS)  VGWuLTS)  II  AMPS)  d(i)  PGEn(r) 

5dJ  * 

5'Jj  * INTEGRATE  CURRENT  aRG  PGmER  pRGM  T=0 . , 1=0. 

60u  200  CQNTIhUE 

6 1 0 CAi_L  RfCPd  I 

62 J Ip(LOGN.EQ.O)  00  iJ  23u 

6 3 j PR  1 1 4 T * "RKPdl  LuuR,  ISC, I" 

j4vj  PRii^i  2 , LUoN , I SC , I 

65o  vjJ  i G 120 

6ju  230  CONTINUE 

o/o  IP  ( (T.LT.TCJT) .Ok . (LCO l . EQ . I ) ) OU  TO  23o 
Odo  EPr  = hLGAl)/bCER 

65 o Ed rLO= ( ( Lm+PML*8S*i ) * 1**2+ ( 1 . - dSR )*  PML*SAiCuR**2 ) /2 . 

/Jo  PR  I NT  * "nGER , Eb  PLU , ELGAl),  E pp" 

/I  o PRInT  d,  EGER, EbP"Li),ELOAJJ,EPP 

72o  LCui'=  1 

7 Ju  235  CONTINUE 

/4o  iP(i-fPG)  250,240,240 

240  CONTINUE  * 

760  IP'(J.PR.EQ.O)  uO  TO  245 
7/u  PRINT  I , T,VC,I,B,L)(2),L>13) 

730  245  CONTINUE 

79u  TPO=TPO+L)TPO 

8 0 0 Ip'(T.GE.TMAX)  GO  TO  200 

dlu  250  CONTINUE 

820  CALL  RKPB2 

830  IF(LDGN.EQ.O)  GG  TO  200 

840  PRINT*"RKPd2  LUGN,  ISC,  I " 

850  PRINT  2 , LDGN , I SC , I 
8 )0  GO  TO  120 
d /O  * 

830  * INTEGRATION  COMPLETE  iO  T=TmaX 
850  260  PRINT* "ELGAG,EpF" 

900  Epp=ELGAG/EGEg 
910  PRInT  3,  ELGAL),  Ep"  p 
920  PRINT*"  " 

930  oG  TO  120 

940  I pURMAT ( I h , dPPd.3,  I PoEl  I .3) 

9dU  2 rGRMATdh  , 1 4,  I PoE  I I . 3 ) 

960  3 pGRMaT ( I H ,IP0EII.3) 

9/0  ENu 


PLu  C aa’I  IS  ) 
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9 J J * 

99u  * CrtLCdLAi  I Qi'4  Uw  JEW  1 \nl  i V ES  ruA  KnHb  rtOuliWc 
I jjJ  SuoHOUriwt  DEW  I \ 
iJIJ  RdnL  I , L,  LL , i-M , idM 
I )eL  J * , I SC , AUu 

IjjJ  COMMOw  i , uT  , v ( 3 ) , u(  J)  ,L,  A 

I J4  J ;,oO,C,w,  IhU.LM,  AM 

IJoJ  J. , d , W , TC  u T , LL , W L 

I Jo J 3. « lUGW , I SC t vb( BMAX , LdCK 

1 j/J  i.MU'J,  JWPM,  HCWM ,DMAC , CUwST,  BDtlLL,  LF’iS 

I -/d  J 4,tiSw , J^L.dWMAX  , AMEt-r 

I _/9J  CaLL  DEW 

I I UO  I rlLJJw  ..(h.  J)  RETURw 

I I I J VTJT* VC- vE 

1120  Ir  ( VTOi'.Li  .0. ) vrJT=D. 

I1SJ  100  D(  1 ) = ( /fOT-A*V(  I ) )/IL+BSh*AML) 

I i 40  D(2)  = VfOT*v<  I ) 

1 I jJ  D(J)  = (LL*DC  1 )+ALxV(  1 ) )*U  I ) 

1 I oJ  * WOT E D(2>  ADD  u(3)  nnt  GEw  aQuER  AWL)  HQnEW  i'G  LOAD  RESdEC i 1 v ELV 
) I / J REi'URw 
1 1 dO  EuD 
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(. 


4dO  * 

44  o * SObROU I h.t:  oEN  CALCULATES  a aS  a ruNCTIOn  D r J AnD  7 AS 

4o0  * A HUNCTION  Of  i AhL)  b 

46o  SUBROUTINE  GEn 

4/o  WEAL  i'I  , 1 , I lsM  ,lM.  t_,  LL 

4dJ  o>,  iscr,  I IT 

49o  Ck,  IC , I J,  ISC,  li 

53J  <j>,  MUO 

5 1 0 COMMON  F , u i , V t d ) , u ( 3 ) , L,  R 

52  o <s,  dO ,C ,N , I dM , LM , lKM 
bdu  Ck,  d,  Vv,  TCUT,  LL,RL 

d4o  tx,  Loun , I SC  , vE , tS.AnX  , LdCn 

bdu  e» , MUO,  dRPM , HCrM , LIMA  G , CUuST , SOriLL , LPFS 

Sou  Ck,  ttbM,  P'AL,  bPMAX,  PMEEr 

b/o  &,SaTCUR 

bdo  DIMENSION  bi (I 0) , 1SCT( lo) , I n ( 10),  70CT(  1 o> , Vj.T(  10) 
i>9o  &, Vi  Si ( 1 O) 
oJo  oAirt  LPRM/I/ 

61  u UAFa  NiAd/// 

0 20  OATa  d I/O . , . o, 1 . , 1 . 9 ,2 . , J . , 4 . , 3*0 ./ 

63  o OAT  a ISCT/O . , 2b.  IEd,4d.lEd,6b.E3,60.Ed,  /o.2Ed,  7o  . Ed , d*0. / 
64j  Oaia  I FT/O .,  2b . I Ld , 4 d . I Ed , 65 • lj , Od .8E3 , 4 7 • cd , d /. 3Ed , 3*0  . / 
6dO  DATa  73CF/U . , /o . , I 34 . , 240 . , 249 . , 370 . , 496 . , d*0. / 

660  OA'iA  V l'f/4*0.  ,24  .,  1 3d. , 225. , d*0 ./ 
o/o  UaTa  Vi  ST/ 4*0. , 5 ., 94 ., 1 bO. , 3*0 ./ 

630  EQUIVALENCE  (1,7(1  )) 

690  IE(LPRM.EQ.O)  GO  TO  90 
/ )o  * 

/lo  * CHANNEL  IDEh i I rICAT lun 

/2o  PR  IwF*  "I  /-400, 2*  1 CHANNEL, QR=0  . " 

730  LRkM=0 

74j  RETURN 

7oo  93  CONTINUE 

760  * , 

7/o  * roLIX  OEuSITY  OuE  Fu  COIL  AUl)  PERMAhEhT  MaGnET 
780  73=0. 

/ 9 0 dS  ti=  1 . 

30o  SAi'CUR=0. 

810  BCL=  1/  IdM 
820  HC=dCL/NUo 

d30  dPM=80HLL*  ( dRPi»i+OMAG*riC*MUO)  /(  UtoAw+dOHLL  ) 

84o  I r ( d PM .LT. dPMAX ) TJ  92 

8d0  I E ( L NTS . NE  • 0 ) GO  fO  91 

doo  PR  Iwi't  "PSaI  " 

d / 0 PR  I Nl  2 , T 

880  91  dP’‘4=3PMAX 

890  SAiCUR=I 

900  dSn=0. 

910  LPTS= 1 


r 


92  o 
940 

940 
9 :>0 
9ju 
9/u 
0)0 
99  j 
I JO0 
191) 
1 J 2 ) 
1 )4  0 
I _/40 
1 oOO 
I jO) 
I J / J 

1 jO) 

I ,9  ) 
1 I uO 
1 1 1 J 
1 I 40 
1 1 J J 
1 1 40 
110) 
1 I oO 
1 1 /0 
I l dJ 
1 1 90 
12o) 
1210 
I 240 
1240 
1 240 
1 2o0 
1 2o) 
12  19 
1 200 

12  90 
I 3o0 

13  10 
1 320 


* r'lV.E  rr  ACCOUNTS  i-QR  rHiooE  t-LuX  ( /vSSUMb  UUl  U t-  CHAN  dUT  LInKS  CUIL3) 
92  4=dCL+bPM*PMErr 
Ir ld.Lf.dMAX ) OL)  TU  9b 
o=BMaX 

Ir (LbCK. EU  . 1 ) oj  iJ  9o 
rR  I if  I * "BMaX.T  ",bllHX,i 

L3Ci\  = I 
vd  CONTINUE 

Hr 


* -jETERATUR  9-1  CdARAChR 1ST IC  AT  CALCULATED  Vy\LUE  Ob  b 

I ER=0 

I oC=  Til  f I (b,NfAd,dT,  I SCI  , 2 , I EH  ) 

I r ( 1 . OE . I CUT ) REiURn 
I f=To f 1 ( b , NTAd , df , I ff . 2 . I ER ) 

VJC  = f.Jri  ( d.NlAd.Bf,  Vl)CT,2,  I EH) 

9 i = ToTl  (B.iirrtb.bT,  Vii  ,2,  I ER ) 

9 TS= IN  f I ( b, NTAd, bf,  Vl'bT,  2 , 1 EH  ) 

Ir(IER.NE.O)  03  TO  129 
I rC 1 1 .LT.O. ) .OR .( I .of . ISC) ) oU  TO  700 
I r ( I .L  f.i'i  ) oO  TO  1 OO 
VA=0 . 

9B=\,fS 

IC=ISC 

1 o=  I T 

J Li  TO  2 00 
1 oO  V A=V  f 
90=  W3C 
IC=I  T 
lu=0. 

* 

* oEnERai OH  vDLI'aoE 

2 00  vO=  VA  + ( 9d-9A) *(  I C-I)/(  IC-  ID  ) 

RETURN 

7oO  LDDN= I 
REiuHu 

/ 20  PRIWf*"frtT|  LHR  IEH , B" 

PRINT  l.lER.b 
STOP 

1 rORMAf ( IH  ,14,1 PbEI 1.3) 

2 rQR MAT ( 1H  ,1 P6EI 1 .3) 

End 
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C . 5 Diet iona ry 
A I ( 2 ) 

AT  (2) 

AV  (2) 

APOLE 

B 

BCL 

BMAX 

BOHLL 

BPM 

BPMAX 

BRPM 

BSW 

BT(10) 

C 

CONST 

D(3) 

DERIV 

DMAG 

DT 

DTPO 


for  Code  COUP  and  Subrout  ines 

Constants  of  linear  equation  for  ampere  turns  per  tesla. 

Constants  of  linear  equation  for  tau  (GFNMAC). 

Constants  of  linear  equation  for  reciprocal  of  volts  per 

turn  per  tesla. 

Area  of  permanent  magnet  pole  face  (GENMAG) . 

Magnetic  flux  density  (tesla). 

Flux  density  due  to  coil. 

Upper  limit  on  total  magnetic  field. 

Absolute  value  of  slope  of  permanent  magnet  circuit  load 

line  on  a plot  of  B vs.  u H. 

o 

Flux  density  due  to  permanent  magnet 

Saturation  flux  density  in  permanent  magnet. 

Residual  Induction.  Flux  density  in  permanent  magnet  when 
H = 0.  With  HCPM,  defines  linear  demagnetization  curve  for 
PM. 

Switch  for  PM  Inductance.  Unsaturated  = 1.  Saturated  = 0. 

Tabulated  values  of  B for  generator  performance  man. 

C/D,  for  magnet  coil.  C/D^  is  a measure  of  coll  cross- 
sectional  area. 

DMAG/ (MU0*IBM)  constant  for  a given  magnet  (i.e.,  C and  N 
specified) . 

Derivative  array: 1 = Current,  2 = Generator  Power  and  3 = 
Load  Power. 

Derivative  subroutine. 

Slope  of  permanent  magnet  magnetization  curve  DMAG  = BRPM/ 
(Ml'0*HCPM)  . (Note  use  of  absolute  value  for  HCPM  is  con- 
sistent with  positive  slope  for  DMAG.) 

Time  Interval.  DT  is  set  at  FT*TAU. 

Time  Interval  for  printout.  (Input  in  milliseconds.) 


C.-l  l 


f 

EBFLD 

EFF 

EGEN 

ELOAD 

FT 

GEN 

HC 

HCPM 

I 

IBM 

icj 
ID  J 

IF.R 

IsC 

IT 

ISC.TT  ( 10)1 
ITT (10)  J 

L 

LBCK 

LCUT 


I.DGN 


LL 


EM 


Energy  stored  in  B field  at  any  instant. 

Energy  transfer  efficiency,  EFF  = ELOAD/EGEN. 

Energy  output  of  generator  integrated  from  t = 0. 

Energy  delivered  to  load.  Integrated  from  time  zero. 

Fraction  of  time  constant  for  one  step  in  time  integration. 

Subroutine  which  computes  MHD  generator  magnetic  field  and 
V-I  characteristic. 

Magnetic  intensity  due  to  coil.  HC  = BCL/MU0. 

Absolute  value  of  intrinsic  coercive  force  of  permanent 
magnet.  Value  of  H in  PM  for  B = 0.  See  BRPM. 

Current  in  generator,  coil  and  load. 

Amperes  per  tesla  for  specified  magnet  coil  (C  and  N 

specified) . 

Characteristic  values  of  current  in  the  MHD  generator.  If 
supersonic  flow  IC  = IT,  ID  = 0;  If  subsonic  flow  IC  = ISC, 

ID  = IT. 

Error  indicator  for  interpolation  routine  TNT1. 

Short  circuit  current  (ISC  is  a function  of  B) . 

Current  in  generator  at  transition  from  supersonic  to  subsonic 
flow  (IT  is  a function  of  B). 

Tabulated  values  of  ISC  and  IT  as  a function  of  BT  which 
characterize  the  MHD  generator.  See  also  VOCT,  VTT,  VTS. 

Inductance  of  coil  plus  load.L  = LM+LL.  Note  LPM  is  treated 
separately. 

Indicator  to  limit  B to  maximum  value  for  which  generator 
characteristics  are  tabulated.  If  B < BMAX,  LBCK  = 0. 

If  B > BMAX,  B is  set  to  BMAX  and  LBCK  = 1. 

Program  switch  to  change  logic  at  T = TCUT.  If  T < TCUT, 

LCUT  =0.  If  T > TCUT,  LCUT  = I. 

Program  indicator  for  error  diagnostic.  LDGN  ^ 0 causes 
print  of  warning  message. 

Inductance  of  load. 

Inductance  of  magnet  coil. 


C-12 


ASSESS 


Program  switch  for  optional  printout  of  variables  as  a 
function  of  time.  If  LPR  = 0,  only  final  results  are  printed. 


Switch  in  subroutine  (JEN  to  print  description  of  MHD  generator 
prior  to  start  of  calculations. 


Switch  in  subroutine  (JEN  to  print  time  at  which  permanent 
magnet  saturates. 


Subroutine  to  compute  magnet  characteristics  for  specified 
values  of  C and  N for  the  coil  plus  permanent  magnet  size 
and  properties. 


Physical  constant.  12.5664E-7  henry/m. 


Number  of  turns  in  magnet  coil. 


Number  of  variables  to  be  integrated  by  subroutine  RKPB. 


Number  of  tabulated  values  for  interpolation  routine  TNT1. 


Inductance  due  to  permanent  magnet  flux. 


Permanent  magnet  effectivenss  factor  to  account  for  fringe 
flux.  Fringe  flux  is  assumed  outside  of  channel  but  still 
linking  the  coil. 


Total  resistance  R = RM+RL. 


Calls  to  library  subroutine  RKPB  for  4th  order  Runge-Kutta 
integration  of  a set  of  first  order  differential  equations. 


Resistance  of  external  load. 


Resistance  of  magnet  coil. 


Coil  Current  at  which  permanent  magnet  saturates. 


Time.  T = 0 at  start  of  MHD  flow. 


Time  constant  of  magnet  coil  alone. 


Time  constant  of  magnet  coil  including  the  effect  of  flux 
due  to  the  permanent  magnet. 


Time  of  cutoff  of  MHD  power.  Current  continues  due  to  circuit 
inductance  (input  in  milliseconds). 


Storage  array  required  by  RKPB. 


Maximum  time  for  calculation  (input  in  milliseconds). 


TNT  1 


TPO 


V(3) 


VA 

VB 


} 


VBM 


VE 

VG 


VOC 

VT 

VTS 


VOCT (10)  ) 
VTT(IO)  ( 
VTS (10)  ' 


Library  subroutine  for  intemolat  ion 

Time  for  which  results  are  printed.  After  printout, 

TPO  = TPO+DTPO. 

Variable  array  1 = I,  2 = Generator  energy  output,  3 = Energy 
delivered  to  load. 

Characteristic  values  of  voltage  on  the  .MUD  generator.  If 
supersonic  flow  VA  = VT ; VB  = VOC.  If  subsonic  flow  VA  = 0, 

VB  = VTS. 

Coil  voltage  per  tesla  for  specified  magnet  coil  (C  and  N 
specif ied) . 

Electrode  voltage  drop. 

Generator  output  voltage.  VG  is  determined  by  generator 
current  and  magnetic  field. 

Open  Circuit  Voltage. 

Generator  voltage  for  supersonic  flow  and  current  corresponding 
to  supersonic  to  subsonic  transition.  See  VTS. 

Generator  voltage  for  subsonic  flow  and  current  corresponding 
to  supersonic  to  subsonic  transition.  The  effect  of  transition 
is  postulated  to  be  a step  decrease  in  voltage  from  VT  to  VTS 
with  no  change  in  current. 

Tabulated  values  of  VOC,  VT  and  VTS  as  a function  of  BT  which 
characteristize  the  MHD  generator.  See  also  ISCTT,  ITT. 
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